Spectroscopic factors of resonance states with the Gamow shell model 
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We provide an investigation of the spectroscopic factors of resonance states in A — 5 — 8 nu- 
clei using two distinct nuclear models, the Gamow shell model (GSM) and the no-core shell model 
(NCSM). Configuration mixing is treated exactly in both GSM and NCSM. GSM employs the 
complex-energy Berggren ensemble, which treats bound, unbound resonance, and non-resonant con- 
tinuum single-particle states on an equal footing. As a result, continuum coupling is taken into 
account in GSM, whereas it is absent in NCSM, where the harmonic oscillator basis is used. We 
first calculate low-lying states of helium isotopes and isotones. The results indicate that GSM can 
more accurately reproduce the low-lying resonance states than NCSM. The spectroscopic factors 
of the resonance states are also computed. Results from the GSM and NCSM are compared and 
show that continuum coupling plays an important role in the calculations of spectroscopic factors 
of resonance states. Overlap functions related to the resonance states are also systematically calcu- 
lated by the GSM and NCSM. It is then demonstrated that wave function asymptotes can only be 


reproduced in GSM, where continuum coupling is taken into account. 


I. INTRODUCTION. 


Large scientific facilities have been upgraded or are 
being constructed, such as the facility for rare isotope 
beams (FRIB), RI Beam Factory (RIBF), the Heavy Ion 
Research Facility in Lanzhou (HIRFL), and High Inten- 
sity Heavy-ion Accelerator Facility (HIAF), which have 
advanced the study of exotic nuclei away from the val- 
ley of stability and propelled it to the forefront of nu- 
clear physics research. The peculiar phenomena that 
distinguish unstable nuclei from those close to the val- 
ley of stability include shell evolution [1, 2], halo and 
skin nuclei [3, 4], exotic decay such as two-proton or two- 
neutron radioactivity [5-8], etc. Such phenomena occur 
as a consequence of valence nucleons primarily occupy- 
ing shells near particle-emission threshold. The study of 
the structure of exotic nuclei has generated great inter- 
est both theoretically and experimentally. Many works 
have already been done on these subjects, particularly 
about halo nuclei [9-12]. However, some questions re- 
main unanswered regarding the coupling of resonance 
and scattering states above particle-emission threshold 
[13], while relevant experimental data are lacking or even 
unavailable. Resonance nuclei located far from the valley 
of stability and close to, or beyond drip-line, form a basic 
laboratory to study the coupling of discrete states with 
scattering continuum states [14]. Traditional shell mod- 
els have difficulties in describing the properties of these 
resonance nuclei. It necessitates the introduction of new 
and powerful models that account for the complicated 
interplay among the different degrees of freedom in those 
weakly-bound and resonance states [15, 16]. 

A typicaly feature of many-body resonance states is 
the strong coupling to the continuum. The coupling be- 
tween the resonance states and the continuum poses the 
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greatest challenges in the theoretical description of un- 
bound nuclei. Various studies have demonstrated that 
continuum coupling has a significant impact on the prop- 
erties of the nucleus in the drip-line regions, such as 
energy spectra [17, 18], asymptotic normalization coef- 
ficients (ANC) [19], spectroscopic factors (SF) [20-22], 
etc. 


In nuclear structure studies, the shell model embed- 
ded in the continuum (SMEC) [16, 23], which treats the 
bound states and the scattering states on an equal foot- 
ing, has been developed these recent years. However, 
SMEC typically only allows one particle to be present in 
the continuum. While it is in principle possible to include 
two particles in the continuum in SMEC, as in the con- 
text of two-proton radioactivity [24, 25], it can only be 
performed with approximations, such as cluster and se- 
quential emissions. In contrast, the Gamow shell model 
(GSM) used in the present work effectively avoids this 
defect and allows multiple particles to be present in the 
non-resonant continuum. GSM allows to perform many- 
body computations using the Gamow-Berggren basis in 
rigged Hilbert spaces [26-29]. The Berggren basis ex- 
tends the Schródinger equation to the complex momen- 
tum space, whereby the scattering contours associated to 
resonance states are discretized with the Gauss-Legendre 
rule. Physical quantities are complex, so that the width 
of resonance states is obtained from the imaginary part 
of energy. A variety of phenomena induced by the prox- 
imity of the nucleon-emission threshold shows that the 
asymptotic part of many-body wave functions must be 
precisely reproduced, which is the case in GSM. GSM has 
thus been used to investigate unbound states in drip-line 
nuclei, such as the conditions of the existence of a 4n 
resonance system [30], the prediction of the ground-state 
energy of 780 [31], and the calculations of energies and 
widths associated to the one-proton and two-proton de- 
cays of !*Mg [18, 32] and to the possibility of a narrow 
resonance ground state in "H [33]. 


SF and its associated overlap function are useful tools 
for the study of the structure of many-body wave func- 
tions of unbound nuclei. SF represents the amount of 
occupancy of a specific single-particle (s.p.) orbit and 
associated shell model configuration in the many-body 
nuclear wave function. It is also a physical quantity 
that connects the nuclear structure to the nuclear re- 
action [34-36]. However, at present, theoretical SFs dif- 
fer significantly from experimental data [36-39]. Stud- 
ies have also demonstrated that continuum coupling pro- 
vide large contributions to SFs in resonance and weakly- 
bound nuclei [22, 39]. As both continuum effects and 
inter-nucleon correlations are taken into account simul- 
taneously in the GSM [40], the asymptotes of many-body 
systems are precisely evaluated. We will employ GSM to 
calculate overlap functions and SFs of resonance states, 
as well as to pinpoint their internal physical mechanisms. 

Ab initio no-core shell model (NCSM) is also consid- 
ered for comparison, in which all nucleons are active and 
interact via a realistic interaction [41]. The basis of har- 
monic oscillator (HO) states is indeed suitable only for 
closed quantum system (CQS) [41], in which only bound 
valence shells are mainly occupied in the s.p. basis. 

After a brief introduction to the theoretical frameworks 
of GSM and NCSM, we compare the results of calculated 
energy spectra and resonance widths of neutron-rich he- 
lium isotopes and their proton-rich isotones with exper- 
imental data. SFs and overlap functions of resonance 
states will be compared using both GSM and NCSM. 
Results allow to investigate the effects of continuum cou- 
pling in resonance states. A summary is provided after- 
wards. 


II. METHOD. 


'The Gamow state, also known as complex-energy res- 
onance state, was firstly proposed by George Gamow in 
1932 to describe o decay [42]. By describing a nuclear 
state in the complex-energy plane, the theory of Gamow 
resonance state provides a natural definition of resonance 
state [42]. Resonance states lie above decay thresholds, 
so that they can decay via the emission of a nucleon, or 
of clusters such as deuteron, ?He/triton, o particles, etc. 

For a given GSM eigenstate, the eigenenergy of the 
many-body state can be written as [43], 


E-E-ir/2, (1) 


where the real part E corresponds to the position of the 
resonance, while the imaginary part I is the resonance 
width, which is related to the lifetime of the state. 

The Berggren completeness relation allows to include 
resonance states in a complete set of one-body states. 
It had been firstly formulated by Tore Berggren and its 
expression is as follows [26] : 


So un(r)un(r’) + I. u(k, r)u(k,r')dk = 6(r — r^), (2) 


where un(r) is a bound or resonance state and u(k,r) 
is a scattering state belonging to the L4, contour in the 
complex-momentum plane. The Berggren basis is indeed 
the extension of the real-energy Newton completeness re- 
lation to the complex plane [44]. Contrary to the New- 
ton completeness relation, which only expands localized 
bound states, the Berggren basis can expand unbound 
resonance states. Within the Berggren representation, 
continuous states are located on the contour L, of the 
fourth quadrant of the complex plane, and the resonance 
states lying between the contour L} and the real-energy 
axis are included in the discrete sum of Eq. (2) along 
with bound states. Continuous states on the contour L4 
are discretized using the Gauss-Legendre method in or- 
der for the completeness relation of Eq. (2) to be used in 
numerical applications. 

Due the presence of complex-energy states in the 
Berggren basis of Eq. (2), the Hamiltonian of GSM is rep- 
resented therein by a complex symmetric matrix, which is 
diagonalized using the complex extension of the Jacobi- 
Davidson method [29]. Resonance states are identified 
in the complex-energy eigenspectrum of unbound states 
formed by both resonant and scattering eigenstates by 
using the overlap method [27, 29]. 

In the ab-initio NCSM, all nucleons in the A-body nu- 
clear system are active, so that the NCSM Hamiltonian 
writes : 


4 vios) so, 


where m is the nucleon mass, View is the realistic two- 
body interaction, to which three-, four-, ... body forces 
can be added in principle. We will only use two-body 
interactions in this work. 

In NCSM, the dimensionality of the matrix explodes 
as the number of nucleon increases, so that a truncation 
of the model space is typically required. In our calcula- 
tions, the standard Nmax truncation is used, i.e. the HO 
energy of basis Slater determinants is at most Nmax hw. 
The convergence of NCSM eigenenergies with respect to 
Nmax can be found in Ref.[45, 46]. As the many-body 
Schródinger equation of Eq. (3) is diagonalized using a 
finite number of HO basis states, whereby Nmax < 20 is 
typical, its eigenstates are deeply bound independently 
of their eigenenergy value. 

The overlap function I;;(r) is defined as: 


Lr) = (Wat t) eei? p, 


J dc 
Mah oj Baca Yun (7); (4) 
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where |W^^) and [uy are the wave functions of the 
A and A — 1 nuclear systems, respectively; / and j are 
respectively the orbital and total angular momenta of the 
considered partial wave in I;(r); al. ej İS the creation op- 
erator associated with the |n;£j) state and is either the 


nj-th state of the HO basis (NCSM) or of the discretized 
Berggren basis (GSM); u,,(r) represents the radial wave 
function of the |n;£j) state. Overlap functions of Eq. (4) 
are independent of the used Berggren basis due to its 
completeness properties, contrary to NCSM, where over- 
lap functions depend on Nmax- 

SF is defined as being the norm of the radial overlap 
function Ip;(r) of Eq. (4): 


Cg = J MOLA (5) 


where C?S is a standard notation for SF. As the GSM 
framework makes use of complex-energy basis states, SFs 
are complex in general, contrary to the real-energy case. 
'The imaginary parts of SFs can be interpreted as the un- 
certainty on real parts generated by the finite life-time 
of the nuclear system [27, 29]. From Eqs. (4,5), the com- 
plex effects of inter-nucleon correlations and continuum 
coupling on the asymptotic behavior of resonance states 
can be directly pointed out [20, 39]. 


III. RESULTS. 


The purpose of this study is multifaceted. First, we 
investigate the energy spectrum and widths of resonance 
states of neutron-rich helium isotopes and isotones with 
GSM and NCSM. Second, we study how SFs and overlap 
functions are affected by threshold effects and coupling to 
the non-resonant continuum. We aim to outline the fun- 
damental features and characteristics of resonance states 
using SFs and overlap functions. 

We will use GSM in the core plus valence particle pic- 
ture. The GSM Hamiltonian will then consist of a one- 
body finite depth Woods-Saxon (WS) potential mimick- 
ing the inert ^He core and of a residual two-body inter- 
action among valence nucleons. In the present work, we 
use for that matter the Gaussian-based Furutani, Hori- 
uchi, and Tamagaki (FHT) residual two-body interac- 
tion, which is composed of central, spin-orbit, and tensor 
parts [47]. The Hamiltonian parameters and model space 
of Ref. [48] are adopted in the present work. The FHT 
interaction has been successfully used in the GSM to de- 
scribe the properties of A — 5 — 10 light neutron-rich 
nuclei [33, 48], the mirror states in !9F/!6N [49], and ra- 
diative capture reactions of A — 6 — 8 nuclei [50]. The 
Coulomb interaction is exactly included in the proton 
space. We will consider in this work weakly bound and 
resonance p-shell nuclei. All partial waves up to l = 3 are 
taken into account. As d3/2 and 7/5, fs/2 partial waves 
bear a large centrifugal barrier, their effect on the asymp- 
totic behavior of the wave function is negligible, so that 
they are expanded with a HO basis. We use the Berggren 
basis for the representation of 51/2, 3/2, 1/2; d5/2 partial 
waves, which form the most important part in nuclear 
wave function asymptotes. The p3/2 and p;;5 partial 
waves are discretized by 60 points along the contour L+, 


TABLE I. The calculated point-proton radius (rp), point- 
neutron radius (rn), and nuclear matter radius (rm) of 59He 
with GSM and NCSM methods, compared with experimental 
data [58-62]. 


6He T» Tn Ty 
NCSM 1.938 2.641 2.429 
GSM 1.868 2.555 2.348 
Exp 1.925(12) [58] 2.45(9) [58] 2.29(6) [58] 
2.47(10) [59] 2.30(7) [59] 
2.51(6) [60] 2.33(4) [60] 
2.61(3) [61] 2.48(3) [61] 
2.66(10) [62] 2.44(7) [62] 
8He f» T qux 
NCSM 1.819 2.570 2.404 
GSM 1.904 2.725 2.545 
Exp 1.807(28) [58] 2.73(9) [58] 2.53(7) [58] 
2.63(9) [59] 2.45(7) [59] 
2.68(5) [60] 2.49(4) [60] 
2.69(10) [62] 2.50(8) [62] 


while the s;/2 and ds/2 partial wave are discretized by 30 
points. The maximal momentum value for the contours is 
Kinase 40 fm !. In performed GSM calculations, three 
nucleons at most can occupy states in the non-resonant 
continuum. These parameters are sufficiently large to 
insure convergence for all the observables of interest. 


In NCSM, we use the Daejeon16 two-body interaction, 
which is based on an SRG-transformed chiral N?LO in- 
teraction [51]. It provides a good description of vari- 
ous observables in light nuclei without three-body forces 
and also generates rapid convergence for energy in ab- 
initio calculations [51]. The HO basis with a frequency 
of fw = 15 MeV is utilized. The model space consists 
of an Nħw space [52], where N < Nmax. We attempt 
to perform computations in large model spaces in order 
to achieve convergent results. However, due to limited 
computing power, maximal truncations are Nmax = 12 
for A = 5 and 6 nuclei, and Nmax = 10 for A = 7 and 8 
nuclei. The results of extrapolating the energy spectrum 
to infinite-size model spaces [46, 53] are shown in Fig. 1. 


The energies of the low-lying states of the ?^5He iso- 
topes and their isotones in ?Li, Be, "B, and 8C nuclei 
are calculated in NCSM and GSM. Results are displayed 
in Fig. 1, along with experimental data [54-57]. As one 
considers mirror nuclei, it is interesting to compare their 
energies and decay widths, which would be otherwise 
identical in the absence of the Coulomb force. Indeed, 
significant differences appear due to the presence of the 
Coulomb Hamiltonian among valence protons in helium 
isotones. A good reproduction of these isotonic states 
is obtained in our GSM calculations. Discrepancy with 
respect to experiment is of a few hundreds of keV when 
using the FHT two-body interaction. On the one hand, 
the calculated ground-state energies of helium isotopes 
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FIG. 1. Low-lying spectra (in MeV) and widths (in keV) 


of the isotopes of ?^^He (upper), and the isotones of ?Li, 
5Be, "B, and °C (lower) The calculations illustrated by 
blue and red lines are provided by NCSM and GSM, respec- 
tively, and the experimental data are taken from Ref. [54-57]. 
Ground-state energies are given with respect to that of “He. 
Shaded areas represent the resonance width in the GSM cal- 
culations and experimental data, with their explicit values 
written above. 


and isotones using NCSM and GSM are both close to 
experimental data. However, the calculated excitation 
energies of excited states in NCSM, on the other hand, 
are generally higher than those of GSM and experimental 
data. The ®SHe isotopes exhibit neutron halos around 
the ^He-like core, which makes a precision measurement 
of the charge radii essential for probing the many-body 
wave function of the neutron halos [63]. We also calcu- 
lated the proton-point, neutron-point, and matter radii 
of 95He with both NCSM and GSM. Results are pre- 
sented in Table. I, along with experimental data. Both 
GSM and NCSM calculations yield good agreement with 
experimental data. 

The GSM decay widths of the neutron-rich helium iso- 
topes agree well with experimental data. As the Coulomb 


barrier is low in light nuclei, unbound states of sizable 
width can form at the proton drip line. GSM accurately 
reproduces the proton-emission widths of the resonance 
states of helium isotones. The difference between calcu- 
lated and experimental widths is indeed smaller than 200 
keV. According to our GSM calculations, ?He(1/2-) is a 
broad (T = 5959 keV) resonance, as well as the "He(5/2-) 
excited state (T = 1831 keV), values which are both con- 
sistent with experimental data, 5570 and 1990 keV for 
?He(1/2-) and "He(5/2-) [55], respectively. 

'The NCSM is diagonalized using a basis of HO states, 
implying that the nucleus is viewed as a CQS. While such 
an assumption is justified for well-bound nuclei, it can no 
longer be applied to unbound states. Weakly bound and 
unbound states bear a low-lying single-nucleon emission 
threshold and a large surface density space dispersion, 
due to a strong coupling to scattering states. It can then 
be expected that GSM will be better at describing reso- 
nance state features than NCSM. 
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FIG. 2. C?S values calculated with NCSM and GSM. Points 
represent the ratio of the GSM results to the NCSM results, 
which is defined as R. Red diamonds correspond to the C?S 
values of neutron, while blue circles are related to the C?S 
values of proton (see right side for specific reaction channels). 
The blue line pinpoints the R = 1 value. SF calculations are 
restricted to the p3/2 and p,/2 partial waves only. 


SF is used as a predictive tool for detailed nuclear 
structure [22, 64, 65], and studies have shown that 
scattering continuum is critical for states near particle- 
emission threshold [22, 66]. To further investigate the 
properties of resonance states, especially the effects of 
continuum coupling, we calculated the overlap functions 
and SFs (C?5 in this paper) of the associated resonance 
states using NCSM and GSM. 

The convergence of C?5S is slower as model spaces in- 
crease in dimension compared to energy in NCSM calcu- 
lations. Added to that, exponential extrapolation of C?.S 
cannot be utilized as for eigenstate energy. Consequently, 
it is necessary to consider the C? and overlap functions 


provided by the largest NCSM model space for compar- 
ison with GSM calculations. In order to investigate the 
differences between GSM and NCSM calculations and to 
demonstrate the continuum coupling effects, the ratio of 
GSM to NCSM results is studied, which is denoted as R 
(Imaginary parts are not considered as they cannot be 
evaluated in NCSM). 

C?S values involving ps /2 and pj? partial waves in 
He, Li, Be, B, and C isotopes bearing A = 6 — 8 nucle- 
ons are calculated and results presented in Fig. 2. Small 
SFs, i.e. for which Re(C?S) < 0.1, are not displayed. 
Calculated SFs are listed in Tabs. II and II. Even 
though SF is a well-defined theoretical concept [67], it is 
model-dependent [68] and can not be observed directly. 
The experimental extraction of SF requires the use of 
the distorted-wave Born approximation (DWBA) [69] or 
eikonal reaction theory [70]. In transfer reaction, SF 
can be obtained by fitting the angular distributions to 
the data points around resonance peaks in DWBA [69]. 
Whereas, in knockout reactions, SF is derived from the 
ratio of experimental and theoretical cross-sections [70]. 
Experimental data provide a SF for 9He(0f) @vp3/2 > 
7He(3/2; ) of 0.512(18) [71], which is close to our GSM 
and NCSM calculations, giving 0.487 and 0.455, respec- 
tively. 

The results of GSM and NCSM differ significantly 
if the R-value departs from the blue line. One can 
see that R ~ 1 in the SFs associated to Ó?He(07) 
Gvpsj > "He(3/2;), SHe(21) vpz > THe(3/2; ), 
and "He(3/2;) ®vp3/2 — 9He(07). Indeed, the corre- 
lated nucleon states of 975He are bound or narrow res- 
onance states in this case, so that continuum coupling 
is reduced. Moreover, in the vast majority of cases, one 
obtains |R — 1| ~ 0.2 — 0.6, thus pointing out strong 
continuum coupling effects. When the A — 1 or/and A 
nuclear states involved in the SF calculation are broad 
resonances, coupling to the continuum is prominent. For 
example, |R — 1| > 0.2 when the A — 1-nucleus is ?Li 
(g.s.) or the A-nucleus is "He(5/2; ), which are broad 
resonances. One of the largest difference between the 
results of NCSM and GSM calculations is indeed ob- 
tained for the reaction ?Li(3/2;) @mp3/2 — 9Be(0j), 
where R ~ 1.6. Similarly, the R value of the reactions 
$He(21) ®vp3/2 > "He(5/2; ) and "He(5/2; ) ®vp1/2 > 
5He(21) is also not close to 1. However, there are excep- 
tions to this rule, with, for example, the reaction ĉBe(0f) 
@vp3/2 > '"B(3/2;). Even though a resonance width of 
about 700 keV is obtained for the unbound "B(g.s.), while 
9Be(g.s.) is also unbound, one has R ~ 1 in this case. 
One can note that a strong continuum coupling does not 
generate either an increase or a decrease of R, as both 
the occurrences R < 1 and R > 1 occur when broad 
resonances are present (see Fig. 2). Illustrations of the 
strong influence of continuum coupling on spectroscopic 
observables in unbound systems can also be found in Ref. 
[72]. 


TABLE II. C?S values issued from NCSM and GSM calcu- 
lations for p3/2 partial wave. R is the ratio between GSM 
and NCSM C?S values. One uses a NCSM model space 
of (NZl, N,.) = (12,12) for A = 5 and 6 nuclei, and 
(Na, NA.) — (10,10) for A — 7 and 8 nuclei. Parent 
and daughter nuclei are denoted by A — 1 and A, respectively. 


C? S(psj3) 
A-1(J")a-1 A (77)a GSM NCSM R 
"He 3/2; He Of  1.850— 0.747 1.572 1.177 
3/27 2t 1.855 + 10.096 1.490 1.245 
1/27 2t 0.138 — 10.015 0.094 1.468 
He Of ‘He 3/2; 0.487-Fi0.101 0.455 1.069 
at 3/2; 1.865 — 10.446 1.931 0.966 
at 5/2; 0.140 — 10.070 0.208 0.675 
"He 3/2; ?He Of  3.222— 10.664 3.384 0.952 
5/27 2t 2.386 + 10.788 1.724 1.384 
"Li 3/27 Be Of 2.421— 10.353 1.540 1.572 
3/21 2t 1.822 + i0.022 1.461 1.247 
°Be Of 7B 3/2, 0.436 +%0.012 0.450 0.969 
at 2.573 — 10.074 1.906 1.350 
"B 3/27 ?C Of 3.788+40.112 3.340 1.134 


TABLE III. Similar to Tab. II, but for C?S of P1/2 partial 


wave. 


C?S(pij2) 
A-1(J™)a-1 A (J")A GSM NCSM R 
"He 1/2; He Of 0.246 — i0.150 0.104 2.363 
$He Of ‘He 1/2; 0.585 — i0.023 0.712 0.822 
at 5/2, 0.821 — 10.005 0.681 1.206 
"He 3/2; 'He 21 0.996 — 10.127 0.788 1.264 


Continuum coupling modifies both SFs and overlap 
functions involving resonance states. As a consequence, 
we will study the overlap functions that correspond to 
the C?S values of Fig. 2. Overlap functions provide de- 
tails on the nuclear structure both inside the nucleus and 
in the asymptotic zone [19, 65]. They are illustrated in 
Figs. 3, 4, and 5, where both real and imaginary parts 
are shown. We only considered overlap functions involv- 
ing the p3/2 partial wave because the latter is dominant 
in the studied nuclei. As the centrifugal barrier is small 
in the / = 1 partial wave, one can expect the A or A— 1 
nucleus, and hence the overlap function, to be extended 
in the asymptotic zone with GSM. This indeed the case 
for all the overlap functions obtained in GSM, which al- 
ways lie above NCSM overlap functions in the asymptotic 
region (see Figs. 3, 4, and 5). This phenomenon is even 
more pronounced when the state of the A or A—1 nucleus 
lies beyond particle-emission threshold. Indeed, due to 
the unbound character of the involved wave function, the 
overlap function can present both negative values in its 
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FIG. 3. Calculated the overlap functions of p3/2 partial wave 
in |A = 5) ® p3;2 — |A = 6) systems with NCSM (dashed 
lines) and GSM (solid lines). The real part for NCSM and 
GSM are shown in the upper panel, and the imaginary part 
of GSM is shown in the lower panel. 


real part and large imaginary parts in modulus, whereby 
an oscillation pattern develops (see the overlap functions 
of the reaction ?Li(3/2; ) @vp3/2 + 9Be(2]) in Fig. 3 
and °Be(0t) &ps/; > "B(3/2,), °He(0t) @vP5/2 > 
7He(3/2; ) in Fig. 4). Conversely, the overlap functions 
obtained in the NCSM are large inside the nucleus and 
rapidly decrease to zero outside the nucleus, which is 
caused by the localized character of the used HO basis 
states therein. 

Nuclei of drip-line regions typically exhibit significant 
isospin-symmetry breaking, such as large mirror energy 
difference of mirror states [73, 74] and abnormal isobaric 
multiplet mass equation [74]. GSM overlap function show 
that there is a large isospin-symmetry breaking in in- 
vestigated mirror partners. For instance, the asymp- 
totic behavior of the nuclei in the reaction ?He(3/2; ) 
@vp3/2 — 9He(21) and in its mirror case, as shown in 
Fig. 3, show noticeable differences. Moreover, the val- 
ues of the overlap function associated to the reaction 
?He(3/2;) ®vp3/2 — °He(2}) is smaller than that of 
its mirror inside of the nucleus. This situation reverses 
outside of the nucleus, whereby the overlap function of 
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FIG. 4. Similar to Fig. 3, but for overlap functions of the 
|A = 6) 8 paja > |A = 7) systems. 


the mirror reaction ?Li(3/2; ) ®mp3/2 + 9Be(21) rapidly 
decreases therein and can become negative. The imag- 
inary parts of overlap functions are also larger in the 
mirror reaction of the proton-rich side, due to the larger 
widths present therein in the A or A—1 nucleus. A sim- 
ilar phenomenon of isospin-symmetry breaking can also 
be seen in the GSM calculations illustrated in Figs. 4 
and 5. On the contrary, the difference between mirror 
nuclei is barely discernible in NCSM overlap functions, 
even though the Coulomb interaction is considered ex- 
actly in NCSM (see Figs. 3, 4, and 5). Thus, the influ- 
ence of continuum coupling in the observables involving 
drip-line nuclei is more significant for isospin-symmetry 
breaking than the simple inclusion of the Coulomb force 
in the many-body Hamiltonian. A large number of HO 
states might be able to generate isospin-symmetry break- 
ing in NCSM results, but the dimensionality of associated 
NCSM model spaces would be enormous, far beyond the 
capabilities of current supercomputers, and that due to 
the slow convergence of SFs and overlap functions with 
increasing Nmax. In fact, only GSM can depict contin- 
uum coupling properly via the explicit use of scattering 
states in the Berggren basis. While the effects of the 
non-resonant continuum can be mimicked by a few HO 
basis states in energy spectra, this is not the case for 
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FIG. 5. Similar to Fig. 3, but for overlap functions of the 
|A = 7) 8 paja — |A = 8) systems. 


observables depending on the structure of nuclear wave 
functions in the asymptotic region, such as SFs and over- 
lap functions. 

Only the C?S value of 8He(07) @vp3/2 > "He(3/2; ) 
has been measured experimentally. In Fig. 6, we present 
calculated overlap function and C?S using both NCSM 
and GSM, along with theoretical results from varia- 
tional Monte Carlo (VMC) and Green's function Monte 
Carlo (GFMC) [75], as well as experimental C?S data 
[71, 76, 77]. Results indicate that the overlap functions 
calculated with GSM and GFMC are more extended than 
those from NCSM and VMC, and that the GSM overlap 
function exhibits the resonance character of 7He(3/27) 
state. However, the C?S values obtained from VMC 
and GFMC are slightly larger than those from NCSM 
and GSM, Nevertheless, all calculated C? values are in 
agreement with experimental data. 

Furthermore, we also performed NCSM calculations 
using chiral N?LO potential with A = 2.0 in SRG and 
chiral N?LO potential, and the resulting C?S values of 
°He(0}) @vp3/2 — 'He(3/2;) are close to each other 
due to the localization of the many-body wave function 
in NCSM calculation. GSM calculations have also been 
performed using a Minnesota interaction with different 
strengths. Results indicate that the calculated C? val- 


ues depend on the one-neutron separation energy (Sn) 
of "He(3/2-) state, and that C?S is close to the results 
obtained with the FHT interaction when the correspond- 
ing Sn is close as well. Therefore, due to near-threshold 
effects, C?S is sensitive to Sy. 


0.5 


= —.He(07) Q vp, — ’He(3/2;) (NCSM) 
$He(01) ® vp — ’He(3/2;) (GSM) 


0.4 — - -ÓHe(07) Q vp, > 'He(3/2)) (VMC) 
$He(07) Q vp, — "He(3/2;) (GFMC) 


x 
T 
& 
— 
E 0.3 
- 
o 
g 
E 02 
ce 
[s] 
E 
9 0.1 
o 
0.0 
0 2 4 6 8 10 12 14 16 18 20 
r (fm) 
0.8 
| E 
0.7 — e 
| GSM & =o & 
e e 
06- GFMC Ss 8| a 
| e a ea x 
N a o $ EE is] 
WY 05 -VMC Š 7 
L e 
04 + NCSM | Q z 
j T 
0.3 - < 
L noce occa 
EXP 
0.2 
FIG. 6. Comparison of the overlap function and C?S of 


9He(0f) @vp3/2 — "He(3/2; ) using NCSM and GSM, along 
with VMC and GFMC calculations [75], as well as experimen- 
tal data [71, 76, 77]. 


IV. SUMMARY. 


The study of unbound resonance states is of great 
significance for the comprehension of the properties of 
atomic nuclei at drip-lines. In this paper, we investigated 
neutron-rich nuclei and their proton-rich mirrors at neu- 
tron and proton drip-lines using the GSM and NCSM. 
The importance of continuum coupling could be precisely 
assessed by comparing nuclear energies, widths, as well 
as SFs and overlap functions. 

First, we calculated the energy spectrum of the he- 
lium isotopes and isotones with GSM and NCSM, where 
the decay widths of resonance states are also provided 
from GSM calculations. We utilized a phenomenologi- 


cal Hamiltonian in GSM, in the core plus valence par- 
ticle picture, with a WS potential mimicking the *He 
and the FHT two-body interaction as residual interac- 
tion among valence nucleons. The Daejeonl6 two-body 
realistic interaction, based on an SRG-transformed chi- 
ral N?LO interaction, is used in the NCSM Hamiltonian. 
The Coulomb force is exactly taken into account in both 
NCSM and GSM Hamiltonians. 

GSM fairly describes the low-lying energy spectrum 
and resonance widths of ?-5He isotopes and °Li, Be, 
7B, and *C isotones. Compared to NCSM results, we 
could already see that continuum effects play an essen- 
tial role in the description of weakly bound and resonance 
states at spectrum level. The asymptotic behavior of res- 
onance states could be precisely investigated in nuclear 
many-body wave functions from the calculation of SFs 
and overlap functions. We have indeed demonstrated the 
effect of continuum coupling in SFs and overlap functions 
of drip-line helium isotopes and isotones. Comparisons of 
calculated SFs in GSM and NCSM show that the contin- 
uum coupling has a significant effect on SFs and overlap 
functions. Indeed, continuum coupling induces changes 
in the SFs of GSM and NCSM, as the ratio of SFs in GSM 
compared to those of NCSM can be both smaller or larger 
than 1, where relative differences can be as large as 60%. 
Added to that, the overlap functions obtained in GSM 
are all extended in space, whereas they are localized in 
NCSM. One could also notice the large imaginary parts 
present in the SFs and overlap functions of GSM when 
unbound states are involved. They provide information 
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